function z = bard(x)
%Bard function
%n=3,m=15,x0=[1,1,1]
%fmin1=8.21487*10^(-3)
%fmin2=f8([0.8406;-inf;-inf])=17.4286
%Author:Ji Lin
y=[0.14,0.18,0.22,0.25,0.29,0.32,0.35,0.39,0.37,0.58,0.73,0.96,1.34,2.10,4.39];
u=1:15;
v=16-(1:15);
w=min(u,v);
z=ones(1,15);
for ii=1:15
    z(ii)=y(ii)-(x(1)+u(ii)/(v(ii)*x(2)+w(ii)*x(2)));
end
z=sum(z.^2);